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PACS 75 . 47 . Lx - Magnetic oxides 

PACS 75 . 47 . Gk - Colossal magnetoresistance 

PACS 75 . 60 . Ch - Domain walls and domain structure 

Abstract. - Several double perovskite materials of the form A2BB'06 exhibit high ferromagnetic 
T c , and significant low field magnetoresistance. They are also a candidate source of spin polarized 
electrons. The potential usefulness of these materials is, however, frustrated by mislocation of the 
B and B' ions, which do not organise themselves in the ideal alternating structure. The result is 
a strong dependence of physical properties on preparative conditions, reducing the magnetization 
and destroying the half-metallicity. We provide the first results on the impact of spatially corre- 
lated antisite disorder, as observed experimentally, on the ferromagnetic double perovskites. The 
antisite domains suppress magnetism and half-metallicity, as expected, but lead to a dramatic 
enhancement of the low field magnetoresistance. 
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^ Double perovskite (DP) materials, of the form 
q A 2 BB'0 6 , are of interest on account of their unusual mag- 
Qnetic and transport properties pi. In particular, some 
'"Hdouble perovskites, e.g., Sr 2 FeMo06, show high ferromag- 
^^netic T c ,~ 420K, half-metallic behavior, and large low 
field magnetoresistance [2|[3] . Other double perovskites, 
CN| e.g., Sr2FeW06, are antiferromagnetic insulators [4], while 
^some show unusual dielectric properties [5]. 

Usually one of the ions, B, say, is magnetic and the 
other (B') is non-magnetic. Magnetic ordering in these 
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materials arises from a combination of strong electron-spin 



coupling on the B ion and electron delocalisation on the 
rHB-O-B' network [6jj8]. The magnetic order, however, is 
K^also strongly affected by the local ordering of the B and B' 

X 



ions |9f|16|. A B-O-B arrangement, for example, favours 
antiferromagnetic locking of the two B moments rather 
Hthan parallel alignment. Most real materials have some 



degree of such 'antisite disorder' (ASD) with mislocated B 
and B' ions. The promise of rich functionality in the DP's 
remains unfulfilled due to this inevitable B, B' mislocation. 

Let us summarise the key observations on S^FeMoOe 
(SFMO), the prototype DP, to establish the tasks for a 
theory of antisite disorder in these materials, (i) Nature 
of defects: there is clear evidence now that B-B' mislo- 
cations are not random but spatially correlated [l7 18 



While ASD suppresses long range structural order, elec- 
tron microscopy LL7] and XAFS [18] reveal that a high 



degree of short range order survives, (ii) Magnetism: the 
ASD brings into play a B-B AF coupling, via B-O-B links, 
that makes neighbouring FM regions antiparallel. This re- 
duces the bulk moment from M max ~ Aps expected in 'or- 
dered' SFMO. The T c is not significantly affected by ASD, 
till very large disorder. The magnetic effects of ASD are 
similar in both single crystals and polycrystals [2]|3{[Io]. 
(iii) Transport: resistivity results are widely different be- 
tween single crystals and polycrystals. Single crystals 3 
show residual resistivity p ~ O.lm^cm, and metallic be- 
haviour, dp/dT > 0. The magnetoresistance (MR) is 
weak, < 10% at low temperature at a field of 5 Tesla. Un- 
fortunately, single crystal data with systematic variation 
of ASD is not available. Polycrystalline samples have been 



studied 10- 14 19 for a wide range of ASD but the trans- 



port in these materials is affected also by grain boundary 
(GB) resistance. The residual resistivity in these samples 
[lO] range from ~ 0.5m^cm for low ASD (M/M max ~ 1.0) 
to ~ lOm^cm at high ASD (M/M max ~ 0.5). 'Ordered' 
polycrystals show dp/dT > 0, while the less ordered ones 
10|[Tl] show dp/dT < 0. The MR can be large, - 40% at 
low temperature and 5 Tesla [IT], and seems to be dom- 
inated by grain boundary effects [13] [l9]. Some results 
indicate a decrease [IT] i n MR with increasing ASD, oth- 
ers show an increase [12] . The effects of antisite disorder 
and grain boundaries on MR have not been deconvolved 
yet. 
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In this paper we build in the correlated antisite disorder 
and study its impact on magnetism and intrinsic transport 
in the double perovskites. We ignore the effect of grain 
boundaries, admittedly important in polycrystals, but un- 
cover a metal-insulator transition and large MR that can 
result from antisite disorder itself. 

Our main results, based on large lattices in two dimen- 
sions (2D), are the following, (i) Correlated antisite disor- 
der leads to structural antiphase boundaries (APB) that 
also act as magnetic domain walls (MDW) at low temper- 
ature, (ii) Growing ASD suppresses the saturation mag- 
netisation, leaves the T c mostly unaffected, but leads to a 
metal-insulator transition in the electronic ground state. 
(Hi) The 'insulating' samples have huge low field MR at 
low temperature since domain rotation 'unblocks' conduc- 
tion paths, (iv) The conduction spin polarisation roughly 
follows the core spin magnetisation in its disorder and tem- 
perature dependence. Overall, moderate ASD seems to be 
beneficial, since it enhances the MR without destroying the 
half-metallicity or suppressing the T c . 



Model and method. 



The ASD arises from an an- 



nealing process, and is spatially correlated 17,18 rather 
than random. Generating antisite disordered configura- 
tions that incorporate the detailed B, B' chemistry and 
capture specific preparative conditions is difficult. How- 
ever, one can model the inherent ordering tendency in the 
B, B' system, and frustrate it via poor annealing to mimic 
the situation in the material. Using the structural motifs 
that emerge [20] we can then solve the coupled electronic- 
magnetic problem through a real space technique. 

We use a binary variable rji to encode information about 
atomic positions. We set rji = 1 for B sites and rji = for 
B' sites. In a non disordered DP the rji will alternate 
along each axis. We will consider spatially correlated dis- 
ordered configurations of 77, discussed later, and for any 
specified {77} background the electronic-magnetic model 
has the form: 



H = 



e B J2 Vifhfi* + cb* X^ 1 ~ m \c 
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Here / is the electron operator corresponding to the 
magnetic B site and m is that of the non- magnetic B' 
site. 6b and eb> are onsite energies, at the B and B' 
sites respectively. e# — e#/ is a 'charge transfer' energy. 
H kin is the electron hopping term: -t x J2(i,j) a ViVjfLfj* 
-*2 E<i, J >( 1 - ~ Vj)ml a m ja -t 3 T,(i,j)*(Vi + Vj - 
^ViV^ifia 171 ^^ + h- c )- The t's are all nearest neighbour 
hopping amplitudes, and for simplicity we set t\ = £2 = 
£3 = t here. is the core spin on the site R^, and 
I Si I = 1. J is the Hund's coupling on the B sites, and 
we use J/t ^> 1. When the 'up spin' core levels are fully 



filled, as for Fe in SFMO, the conduction electron is forced 
to be antiparallel to the core spin. We have used J > 
to model this situation. For the present study we have set 
[21] the effective charge transfer energy e# — J/2 — e#/ =0. 

If two magnetic atoms are on neighbouring sites, we 
also have to consider the antiferromagnetic (AF) coupling 
between them. Much of the physics of these materials 
arises from the competition between delocalisation driven 
ferromagnetic exchange and this B-B superexchange. We 
set JafS 2 jt = 0.08, based on the T/v scale in SrFe03. h 
is an applied magnetic field in the z direction 

The model above treats the variables rji as given. What 
governs the distribution of ionic positions {rji}? The ef- 
fective interaction between ions has a 'direct' component 
(Lennard- Jones, say) that we label Vi on . There is also the 
indirect interaction via the electronic-magnetic degrees of 
freedom in H{rji}. The total Hamiltonian, including all 
degrees of freedom, would be: H tot = H{rji} + Vi on {rji}. 
To obtain the effective interaction among the {77^} one 
should trace out the electronic and magnetic variables. So, 
the effective potential V e ff {rj} controlling the positional 
order is V eff {rj} = -(l//3)log j VSiTr {firn} e-^ Htot . Un- 
fortunately there is limited information about V^ on , and 
the trace is computationally demanding. Thankfully, the 
experiments themselves suggest a way ahead. 

Firstly, the structural degrees of freedom 'freeze' at a 
temperature T ord ~ 1000K > T c ~ 400K, where T c is 
the magnetic ordering temperature, so we need not bother 
about the 'feedback' of magnetic and electronic degrees of 
freedom on the {rji} ordering. Secondly, the tendency is 
to order into an alternating pattern, frustrated by short 
annealing. So, as a first attempt, we can try out a simple 
short range model with the same tendency. Concretely, we 
use a 'lattice gas' model V e ff{rj} = —V^2^rji(l — rjj) 
with V > being a measure of the ordering tendency. 
The ground state in this model would correspond to 
rj = 1,0,1,0,.. along each axis, z.e, B, B', B, B',.. This 
approach tries to incorporate the effect of complex inter- 
actions between the A, B, B' and O ions (as also the elec- 
trons) into a single parameter V. 

We explored [20] the consequences of different annealing 
protocols on this model, to examine the consequences of 
imperfect annealing. The lattice gas model has a finite 
temperature transition in both 2D and 3D. We work in 
2D since it allows access to large sizes and is easier to 
visualise. 

One can variously characterise the B-B' structures that 
emerge. We use the following simple indicators, to keep a 
close correspondence with the experimental work, (i) The 
fraction of B (or B') atoms that are on the wrong sub- 
lattice, call this x, and the structural 'order parameter' 
5 = 1 — 2x. (ii) The degree of short range order, charac- 
terised by the probability, p, of having nearest neighbour 
pairs that are B-B'. (Hi) The correlation length £ in these 
structures, computed from the width of the ordering peak. 

For a given {rji} configuration we need to solve for the 
magnetic and electronic properties. Since the background 
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Fig. 2: Colour online: (a). The ferromagnetic peak in the struc- 
ture factor, M 2 , where M is the magnetisation, for different 
degrees of antisite disorder, (b) The 'd.c resistivity', p(T) for 
varying antisite disorder. Results are on a 40 x 40 lattice, av- 
eraged thermally and over 10 copies of disorder. 



Fig. 1: Colour online: We show four families of antisite dis- 
ordered configurations (top to bottom) generated via succes- 
sively poorer annealing of the lattice gas model. We plot 
(rji - l/2)e t7r ( Xi+yi \ For a perfectly ordered structure g(Ri) 
is constant. The patterns along a row are different realisations 
of ASD within each family. The average structural order pa- 
rameter (see text) has values S = 0.98, 0.76, 0.50 and 0.08 as 
we move from top to bottom. Lattice size 40 x 40. 



involves strong disorder, and the electron-spin coupling, 
J, is large, we use an exact diagonalization based Monte 
Carlo (ED-MC) technique. This uses the Metropolis algo- 
rithm where a spin update, —> S • is accepted or rejected 
depending on AE/k B T, where AE = E(S f i ) - E(Si). In 
principle we should diagonalise the full system every time 
an update is attempted. The cost, for a large system, is 
prohibitive, so we employ a method [22] where we diago- 
nalize a cluster Hamiltonian built around the update site. 

Electronic properties are calculated after equilibration 
by diagonalizing the full system. The optical conduc- 
tivity is calculated via the the Kubo formulation. The 
'dc conductivity' is the low frequency average [23], <j^ c = 
(1/Acj) J Q Au; a(uj)duj, where a(u) is thermal and disorder 
averaged, and Auj ~ O.Obt. Our 'dc resistivity' is the in- 
verse of this <Jdc- The spin resolved density of states D a (u) 
is calculated from the single particle Green's function. 



Earlier work. — The theoretical effort till now |15||16| 
has focused on uncorrelated antisite disorder, and quanti- 
fied the impact of such ASD on magnetic properties. One 



of them 15 is based on a variational scheme in large 3D 



systems, and quantifies the doping and antisite concentra- 



tion dependence of magnetic properties. The other 16 
uses a classical spin model for magnetism, and studies the 
critical properties. These calculations set the reference 
for magnetic properties but have the obvious limitation 
that: (i) they use randomly located antisites, (ii) they do 
not clarify the electronic properties, and (Hi) the estimate 



of localisation effects that arise from structural/magnetic 
domains, and the possible MR, remains unexplored. 

Results. — 

ASD configurations. The ASD configurations can be 
readily generated through a Monte Carlo on the lattice 
gas model. Below an ordering temperature T or d ~ 0.7V 
(in 2D) the model exhibits long range B-B' order provided 
one anneals long enough. We quench the system from 
high temperature (random B, B') to T ann < T or ^ but 
deliberately anneal it for a short time, preventing equilib- 
riation. The details have been discussed [20] in an earlier 
paper. The four families we chose, see Fig.l, have a frac- 
tion of mislocated sites: x = 0.01, 0.12, 0.25, 0.46. The 
patterns, however, are strongly correlated. Even in the 
most disordered samples (lowest row in Fig.l) where the 
likelihood of any site being B or B' is ~ 0.5, if a site 
is B, say, there is a high likelihood that its neighbours 
will be B' (and vice versa). Following recent experimental 
work 18 , we calculate the probability p of having B-B' 
nearest neighbours as a measure of short range correla- 
tion. In a perfectly ordered sample this would be 1 (the 
B, B' alternate) while in a completely disordered sample 
this is 0.5. For an uncorrelated B, B' distribution this is 
Puncorr = x 2 + (1 — x) 2 = (1/2)(1 + S 2 ) . If we only knew 
S, these would lead to p U ncorr = 0.98, 0.79, 0.63, 0.50 
for the x that we have used. The values that actu- 
ally emerge from analysing our correlated patterns, are 
p corr ~ 0.98, 0.97, 0.95, 0.86. Even the most disor- 
dered samples have a high degree of short range order. A 
Lorenztian fit to the B-B' structure factor, of the form 
W(q) ~ r7(fec - tt) 2 + (q y - tt) 2 + r 2 ), yields 
f - 6.6, 5.9, 4.8, 3.6. 

Disorder dependence at T = 0. Let us examine the ef- 
fect of the ASD on the magnetic properties. Suppose the 
fraction of mislocated B, B' sites is x, and the structure is 
organised into domains such that the ratio of 'perimeter' 
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Fig. 3: Colour online: Temperature dependence of short range 
magnetic correlations. The left column shows ASD configura- 
tions, one from each disorder family. S = 0.98, 0.76, 0.50 and 
0.08, top to bottom. The 2nd, 3rd and 4th panel along each 
row is a map of a spin overlap factor, gi, in a Monte Carlo 
snapshot, gi — So. Si, where So is the left lower corner spin in 
the lattice. The temperatures are T/t = 0.03, 0.05 0.07. 
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Fig. 4: Colour online: (a)&(c) Field dependence of magneti- 
sation, M 2 , (b)&(d) field dependence of the resistivity, nor- 
malised to h = 0. The left panels are at 'low temperature', 
T/t = 0.03, the right at 'high temperature', T/t = 0.07. 



to 'bulk' sites is small. The AF coupling between adjoin- 
ing domains would polarise them antiparallel, and the net 
moment at T = 0, h = would be proportional to the 
volume difference of up and down domains. We should 
have M(T = 0, h = 0) - (1 - x) - x = 1 - 2x = S so 
M 2 = S 2 . Given our S, these are 0.96, 0.58, 0.25, 0.01 in 
almost perfect correspondence with the T — > values in 
Fig. 2. (a). The elaborate calculation arrives at an obvious 
answer. The onset temperature for magnetic order seems 
to be insensitive to the ASD, i.e., the intra-domain order 
sets in at T ~ the bulk T c . While our answer for the sup- 
pression of magnetisation is M ~ 1 — 2x, a 3D calculation, 
with uncorrelated disorder, had found [15] M ~ 1 — 1.9x. 

The transport in these background is shown in Fig.2.(b). 
From weak to intermediate ASD the temperature depen- 
dence of p(T) remains similar, with a sharp drop near 
T c . The only effect of increasing ASD is an increase in 
residual resistivity. It is as if there is a temperature inde- 
pendent 'structural' scattering that gets added to the T 
dependent magnetic scattering. At large ASD, however, 
this correspondence breaks down, the T = 'resistivity' 
is very large (and grows with growing system size) and 
dp/dT < 0. There seems to be a metal-insulator transi- 
tion between S — 0.50 and S = 0.08 

To create an understanding of this let us focus on T = 0, 
where the magnetic configuration is simple (collinear). 
The 4- s P m electrons inhabit the t core s P m domains, 
and vice versa. The conductance arises from the inter- 
penetrating parallel channels for up and down spin elec- 



trons. One could call it "complementary percolation" . Let 
us identify up electrons with the 'majority' phase, and 
down with the 'minority' domains. The net conductivity 
°~tot(S) = (Tmaj(S) + Cmin(S). While this reduces mono- 
tonically with reducing S in our data, Fig.2.(b), the con- 
ductivity also depends on p CO rr- Even at S = 0, one could 
increase a tot systematically by increasing p corri i.e, re- 
ducing the fragmentation of the conduction paths. Weak 
localisation effects, etc, in two dimensions could show up 
at much longer lengthscales, but at a given cutoff size the 
trend above would survive. 

Temperature dependence:. The ASD configuration is 
T independent so the primary sources of T dependence on 
transport are (i) the weakening of AF locking across the 
domain boundaries, and (ii) fluctuations about the FM 
state within a domain. The first effect enhances the con- 
ductivity, while the second serves as a source of scattering. 
Their relative importance depends on a(T = 0). For weak 
disorder (large S) one is far from the percolation threshold 
and the decrease in a due to intra-domain magnetic scat- 
tering is larger than the enhancement from inter-domain 
tunneling. However, by the time S = 0.50 there is already 
a weak upturn in p as T — >• 0, the intra-domain effect is 
visible, and this becomes the dominant effect as S — >• 0. 

An analysis of the spatial spin-spin correlations illus- 
trates the AF locking of domains at low temperature and 
how this weakens with increasing T. The first column 
in Fig. 3 reproduces one set of ASD configurations from 
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Fig. 5: Colour online: Field dependence of magnetic spatial 
correlations. We show the usual gi = So. Si, defined earlier. 
The left column shows the ASD domains, the central column 
shows gi at h/t = 0, and the right column is for h/t = 0.001. 
The temperature is T/t = 0.03. 

Fig.l (the first column). The next three columns show 
the the magnetic overlap gi = So- Si, where So is the 
lower left corner spin in each configuration, for MC snap- 
shots at T/t = 0.03, 0.05, 0.07. These pictures would 
correspond to 'magnetic domains' if the patterns survived 
thermal averaging. The low temperature snapshots corre- 
spond closely to the ASD pattern. The antiphase bound- 
ary (APB) and the magnetic domain wall (MDW) pattern 
coincide at T/t = 0.03. At T/t = 0.05, however, close 
to the bulk T c , there is no correlation between the APB 
and the gi pattern. There is significant core spin overlap 
across the boundary, and large fluctuation, overall, in spin 
orientation. This bears out the transport mechanism we 
suggested in the preceding paragraph. 

Field dependence. The field dependence of magneti- 
sation and resistivity is shown in Fig. 4., at relatively low 
temperature, T/t = 0.03, in (a)-(b), and high tempera- 
ture, T/t = 0.07, in (c)-(d). Three energies play out when 
ft ^ 0: (i) the bulk Zeeman cost of the 'minority' domains, 
~ hVmin, where V m in is the volume of the minority phase, 

(ii) the interfacial AF energy, ~ JafV(1 — p CO rr): where 
1 — Pcorr is the fraction of AF bonds on the lattice, and 

(iii) the gain in electronic kinetic energy on removal (or 
rotation) of MDW's. (i) and (iii) prefer domain alignment 
while (ii) prefers to retain domain walls. In a 'spin only' 
model (iii) would be absent. This delocalisation energy 
gain serves to reduce the field at which domain rotation 
can occur. 




Fig. 6: Colour online: Half-metallicity, estimated as (D^(cf) — 
D±(ef)) I '(D^(cf) + Di(cf)), for varying temperature and dif- 
ferent degrees of antisite disorder. 

At low T, Fig. 4. (a)- (b), the ordered samples have a high 
degree of magnetic order, so the field induced increase in 
M, and the decrease Ap/p(0), where p(0) — p(h = 0) and 
Ap = p(0) — p(ft), is quite small. The low T low field 
response is, however, dramatic for low S samples. These 
samples have M(h = 0) ~ 0, and a large p(0) due to the 
fragmented (spin selective) conduction path. A field as 
small as h/t ~ 0.001 leads to M 2 - 0.1, so M ~ 0.3. The 
corresponding impact on spin correlations is shown in the 
lowest row in Fig. 5, where the MDW pattern is strongly 
affected by ft. While the domain rotation effect is visible 
both for S = 0.50 and S = 0.08, the less disordered sample 
had a larger conductivity at ft = so the fractional change 
is much larger for S = 0.08. 

At high temperature, Fig.4.(c)-(d), the domains cease 
to exist and conductance gain from domain rotation is 
irrelevant. In the large S samples there are few AF links so 
the applied field just suppresses the magnetic fluctuations 
leading to large Ap/p(0). In the most disordered samples 
there are (1 —p CO rr)/^ °° 7% of AF bonds. Although there 
are no domains, these act as a source of scattering. The 
gain in conductivity is slower in the disordered samples 
compared to the more ordered ones. 

Half metallicity. These systems are unusual because 
at T = within each domain the conduction electron 
has only one spin polarisation but averaged over the sys- 
tem both t and 4- electrons have density of states at ep- 
A local probe, with probe area < £ 2 , will allow only 
spin polarised tunneling, while a probe averaging over 
domains will see both D^{ep) and D^{ep). Fig. 6 shows 
(D^(cf)— Di(cf)) / (D^(€f)+Di(€f)) as a measure of half- 
metallicity. It is unity only in the absence of ASD and at 
T = 0, and in general has a behaviour that broadly mimics 
the behaviour of the core spin magnetisation, Fig. 2. (a). 

Discussion. — There are two issues we want to touch 
upon, to relate our work to real double perovskites. 
(a). The role of dimensionality: it is well known that lo- 
calisation effects are stronger in 2D compared to 3D, so we 
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ran this entire calculation on a 12 3 system to check out the 
trends in transport. In particular we confirmed that there 
indeed is a sharp rise in the residual resistivity (although 
possibly no insulating phase) with increasing ASD. The 
low temperature upturn in p(T) is present, but weaker, 
even in 3D. (b). Role of grain boundaries: in the absence 
of a chemical characterisation of the grain boundary ma- 
terial, and an electronic model for the GB, it is hard to 
construct a comprehensive theory. However, since grain 
size, Iq » £, it should be possible to study the role of 
APB's and MDW's via non contact probes that focus on 
a single grain. 

Conclusion. — We have studied a double perovskite 
model on antisite disordered backgrounds with a high de- 
gree of short range correlation. In this situation, the an- 
tiphase boundaries coincide with the T = magnetic do- 
main walls. Growing ASD reduces the low field magne- 
tization, destroys the half-met allicity, and leads to a low 
temperature metal-insulator transition. While these are 
disadvantages, we also note that the ferromagnetic T c is 
only weakly affected by moderate ASD and the low field 
magnetoresistance is dramatically enhanced by disorder. 
Our real space results allow an interpretation of these in 
terms of the domain pattern, the effective exchange, and 
the short range magnetic correlations. They are also con- 
sistent with explicit spatial imagery from recent experi- 
ments. The 'ultra-grain' effects highlighted here would be 
directly relevant to single crystals, and define the starting 
point for a transport theory of the polycrystalline double 
perovskites. 
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